O and ^^V NMR for the zigzag spin-1 chain compound CaV204 
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^^V NMR studies on CaV204 single crystals and ^^O NMR studies on ^^O-enriched powder 
samples are reported. The temperature dependences of the ^^O NMR line width and nuclear spin- 
lattice relaxation rate give strong evidence for a long-range antiferromagnetic transition at Tn = 78 K 
in the powder. Magnetic susceptibility measurements show that Tn = 69 K in the crystals. A zero- 
field '^^V NMR signal was observed at low temperatures (/ « 237 MHz at 4.2 K) in the crystals. The 
field swept spectra with the field in different directions suggest the presence of two antiferromagnetic 
substructures. Each substructure is coUinear, with the easy axes of the two substructures separated 
by an angle of 19(1)°, and with their average direction pointing approximately along the 6-axis of 
the crystal structure. The two spin substructures contain equal number of spins. The temperature 
dependence of the ordered moment, measured up to 45 K, shows the presence of an energy gap 
Ea in the antiferromagnetic spin wave excitation spectrum. Antiferromagnetic spin wave theory 
suggests that Ec/kB lies between 64 and 98 K. 

PACS numbers: 76.60.-k, 75.50.Ee, 75.25.-l-z, 75.10.Pq 



I. INTRODUCTION 

Frustrated magnetic systems have attracted a lot of 
research interest because such systems often exhibit 
interesting low temperature properties i The zig-zag 
spin chain with antiferromagnetic interactions between 
nearest- and next-nearest-neighbors is about the most 
simple frustrated system. In a zig-zag spin chain system 
with spin S — 1, the ground state phase diagram (at 
temperature T = 0) as a function of XXZ anisotropy 
and ratio between nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) interactions exhibits six differ- 
ent phases.- '"^ In addition to two Neel ordered phases and 
two phases with a Haldane gap, there exists a large phase 
region called a gapless chiral phase where the chirality ex- 
hibits long range order without accompanying spin order, 
and a small phase region where there is a gapped chiral 
phase. 

CaV204 is a possible candidate for a zig-zag spin S = 1 
chain system. ^'^ It has an orthorhombic crystal structure 
(space group Pnam) at room temperature as shown in 
Fig. [TJ Vanadium moments at two crystallographically 
inequivalent sites respectively form two inequivalent zig- 
zag spin chains along the c-axis. In one of the two chains, 
the distances between NN and NNN vanadium atoms are 
3.01 and 3.08 A, respectively, while in the other chain, 
these two distances are 3.01 and 3.06 A, respectively. 
The smallest interchain vanadium distances are 3.58 and 
3.62 A,^ Thus one might expect a much smaller inter- 
chain coupling as compared to NN and NNN interactions 
within the chain. 

Previous magnetization and ^^V NMR studies of 
CaV204 suggested that the ground state of the system 
might be a gapless chiral phase. However, this find- 
ing contradicts earlier powder neutron diffraction studies 
which showed that the system is in an antiferromagnetic 




FIG. 1: The orthorhombic crystal structure of CaV204 at 
room temperature. Large spheres - two inequivalent V sites, 
small dark spheres - Ca, small light spheres - O. The crystal- 
lographic a and b axes are along the vertical and horizontal 
directions in the plane of the page, respectively. The c axis is 
perpendicular to the page. The cuboid denotes the size of a 
unit cell. 



state at 4.2 K;^ The neutron measurements indicated 
a magnetic unit cell in which the b and c lattice con- 
stants are doubled and the spin directions in each chem- 
ical unit cell are reversed relative to their orientations 
in neighboring chemical unit cells along the b- and c- 
axes. Each chemical unit cell contains 8 vanadium spins 
and the magnetic spin structure within each chemical 
unit cell could not be uniquely determined. The neu- 
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tron diffraction pattern was found to be consistent with 
three different collinear models with spins parallel to the 
6-axis. By assuming the same spin moment at all vana- 
dium sites, the magnetic moment of each vanadium spin 
was determined to be 1.06(6) /ib, only about half the 
value expected for a vanadium spin with y-factor g ~ 2 
and 5' = 1.— The presence of a low temperature antifer- 
romagnetic phase is also supported by recent magnetiza- 
tion measurements on annealed CaV204 single crystals, 
which showed a clear signature of an antiferromagnetic 
phase transition at temperature Tn = 69 K.^ 

In order to obtain further evidence of the magnetic 
phase transition and to study the magnetic properties 
in the ordered state, we performed ^^O and ^^V NMR 
studies on ^^0-enriched powder and single crystal sam- 
ples, respectively. The ^^O NMR spectrum and the 
nuclear spin-lattice relaxation rate measurements give 
strong evidence for a magnetic transition at Tn = 78 K 
in the powder sample. In contrast to early ^^V NMR 
measurements)^!^ we could not detect a ^^V NMR signal 
in the powder sample around the normal Larmor fre- 
quency in an applied field of 1.67 T in the temperature 
range of 4.2 K < T < 296 K. Instead, we observed a 
zero-field ^^V NMR signal at T < 45 K (/ « 237 MHz 
at 4.2 K). 

The zero-field ^^V NMR signal is observed because of 
a strong local field at ^^V nuclear sites (i?ioc = 21.2 T at 
4.2 K) in the ordered state. The local field arises mainly 
from interaction between nuclei and vanadium core elec- 
trons, which are polarized by the ordered 3c? electronic 
spinsjSiii^ This local field points antiparallel to the di- 
rection of the local electronic spin moment. By studying 
how the resonance frequency changes as a function of the 
direction and magnitude of the applied magnetic field, 
one can obtain information on the vanadium spin struc- 
ture, as will be demonstrated below. We measured the 
temperature dependence of the ordered moment to study 
the anisotropy gap of the antiferromagnetic spin wave ex- 
citations. We also attempted to measure the temperature 
dependence of the V nuclear spin-lattice relaxation rate 
l/Ti. However, due to the very broad line and the pres- 
ence of nuclear quadrupole splitting (the nuclear spin of 
^^V is / = 7/2), the relaxation curves are highly nonex- 
ponential and depend strongly on the saturation condi- 
tion. Thus, a reliable measurement of the ^^V l/Ti(T) 
was not achieved. 

The remainder of the paper is organized as follows. 
Experimental details are explained in Sec. [ill ^^O and 
^^V NMR results are presented in Sees. IIIII and IIVI re- 
spectively. In Sec. |Vl we give a summary of the main 
results of the paper. 



II. EXPERIMENTAL DETAILS 

Polycrystalline single phase CaV2 04 (sample an-2- 
116) was synthesized via the solid-state route by reacting 
V2O3 (99.995%, MV Labs) with CaO obtained by calcin- 



ing CaCOg (99.995%, Althaea Chemicals) at 1100 °C. 
The chemicals were ground inside a He glove-box, then 
pressed and sintered at 1200 °C for 96 hours in flowing 
4.5% H2-He with intermediate grindings. Phase purity 
was confirmed by powder X-ray diffraction (XRD) on a 
Rigaku Geigeflex diffractometer using Cu Ka radiation in 
the 29 range of 10°-90°.-& 72.1 atomic% ^'^O-enriched oxy- 
gen (MSD Isotopes) was used for ^^0-enrichment. About 
1 g of the precursor CaV204 was placed in a Pt foil-lined 
alumina boat in an evacuated furnace tube, which was 
then preheated to 750 °C under dynamic vacuum. The 
pumping line was then closed and the tube backfilled 
with the i^O-enriched O2. The mass gain on oxidation 
indicated a nominal composition of CaV205.94. This was 
placed in flowing 4.5% H2-He and reduced as before to 
CaV204 (sample an-2-180E). The final ^^O content of 
the enriched CaV204 was about 25%. Powder XRD was 
used to confirm that the sample was single phase. 

^^O NMR measurements were performed utilizing a 
phase-coherent pulse spectrometer in applied fields of 3.0 
and 4.7 T. The typical 7r/2 pulse length is 6 /is. The echo 
signal was produced by a sequence of a tt/2 and a 7r/3 
pulse, which produce the maximum echo signal intensity. 
The separation between these two echo generating pulses 
was 40 us. The ^'^O NMR spectra were measured by ei- 
ther Fourier transform of half the echo signal or by plot- 
ting the area of the echo as a function of the rf frequency 
(frequency sweep). The nuclear spin-lattice relaxation 
rates were measured by monitoring the recovery of the 
echo intensity following a comb sequence of 7r/2 satura- 
tion pulses. Static magnetization versus temperature was 
measured in a Quantum Design SQUID magnetometer in 
a field of 1 T and in the temperature range 5-100 K to 
confirm the low temperature magnetic behavior and the 
ordering temperature Tn = 78 K. The magnetic suscep- 
tibility of the powdered "'^^0-enriched sample is shown in 
Fig. O The transition temperature is revealed by a small 
kink in the x(T) data at Tn. 

Two CaV204 crystals were used in V NMR measure- 
ments. Crystal ^1 was grown in an optical fioating zone 
furnace while crystal #2 was grown using a tri-arc crys- 
tal pulling method. ^'"'^1 The sizes of the crystals #1 and 
#2 are about 5 x 5 x 10 mm'^ and 1x1x2 mm'^, re- 
spectively. Both crystals were annealed at 1200 °C under 
4.5% H2-He fiow and the antiferromagnetic phase tran- 
sition temperatures of the annealed crystals were found 
from magnetic susceptibility and heat capacity measure- 
ments to be 69 K.— The difference in Tn between the pow- 
der and single crystal samples may be related to different 
unit cell volumes of these samples. The unit cell volumes 
in the powder and single crystal CaV204 samples are 
296.0(5) and 298.0(5) A^, respectively. The reason for 
this difference is, however, unknown at present. A sim- 
ilar correlation between Tn and the lattice volume was 
observed in EuCu2Ge2.— The magnetic susceptibility of 
crystal #2 is shown in Fig. [3] with the field along a and b 
directions. The antiferromagnetic transition temperature 
Tn is clearly seen as a bifurcation in the susceptibilities 
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one of the two magnet pole caps. The difference be- 
tween the measured field and the field at the position 
of the sample (measured by the resonance frequency of 
protons in water) was less than 0.005 T over the whole 
field range. Individually, the crystals were placed inside 
copper solenoid coils with the crystal c-axis parallel to 
the coil axis. Measurements of crystal #1 involved rota- 
tion of the field in the a-b and b-c planes. Measurements 
of crystal #2 involved rotation of the field in the a-b 
plane. The rotation of the field was achieved by rotating 
the cryostat together with the crystal about the crys- 
tallographic axes perpendicular to the field plane. The 
misalignment between the rotation axis and the intended 
crystallographic axis is estimated to be less than 5°. 



FIG. 2: Magnetic susceptibility x versus temperature T of 
the ^^O-enriched CaV2 04 powder sample measured in a field 
of _ff = 1 T. The vertical arrow indicates the position of the 
antiferromagnetic transition temperature Tn = 78 K. 



III. 




H = 1 T 



an-2-50-c5, CaV O 

2 4 

annealed triarc crystal 
(NMR crystal #2) 



a//H 
b//H 



100 200 
Temperature (K) 



300 



FIG. 3: Magnetic susceptibility x versus temperature T of 
CaV2 04 crystal #2 measured with the applied field _ff = 1 T 
along the a and b directions. The vertical arrow indicates the 
antiferromagnetic transition temperature Tn — 69 K. The 
measurements were carried out under either field-cooled or 
zero-field-cooled conditions, as indicated. 



along the two directions. We note that, when the field 
is along the b direction, a splitting between zero-field- 
cooled and field-cooled susceptibility is observed below 
T = 20 K. 

A search for a zero-field ^^V NMR echo signal was per- 
formed at 4.2 K and was found to be located close to a 
frequency of 237 MHz at that temperature. The echo 
was produced by a sequence of two pulses with the same 
pulse length, which was typically 4 /is and about half the 
length of a 7r/2 pulse. The separation between the pulses 
was fixed to 18 /zs. ^^V NMR spectra were measured by 
plotting the echo intensity as a function of the magnetic 
field. A variable magnetic field from to 2.0 T was pro- 
duced by an electromagnet. The value of the magnetic 
field was measured by a Hall magnetometer attached to 



^^O NMR IN POWDER SAMPLE OF 
CaVaOi 



Figure |4] displays the ^''O NMR spectra for the ^''O- 
enriched powder sample of CaV204 in = 3 T at three 
different temperatures. The spectrum at T = 296 K was 
obtained via Fourier transform of half the echo signal 
while the spectra at T = 80 and 77 K were obtained 
by frequency sweep. ^^O nuclei have spin 5/2 and thus 
possess a nonzero electric quadrupole moment. Since the 
local environments of all oxygen sites do not possess cu- 
bic symmetry (point group to), one expects a quadrupole 
splitting of the ^''O resonance frequencies. We attribute 
the lack of a powder pattern of the first order quadrupole 
splitting in the observed spectra to a smaller quadrupole 
splitting compared to the magnetic broadening of the 
spectra. The positions of the spectrum peak position 
show negligible field dependence. The chemical shift K 
of the peak position is given by K = AhfXj where Ahf 
is the isotropic hyperfine coupling constant and x is the 
magnetic susceptibility. Using the data in Figs. [5] and 
Uland additional x{T) data (not shown), we obtain the 
upper limit Am < 1.7 kG//j,B- 

The line gets broader below Tn with decreasing tem- 
perature. As shown in Fig. 21 the absorption line at 
T = 77 K is much broader than the lines at T = 80 and 
296 K. As will be further shown below, this broadening 
is a signature of an antiferromagnetic phase transition at 
Tn = 78 K, where the nuclear spin-lattice relaxation rate 
l/Ti exhibits a peak. As the temperature approaches 
the phase transition temperature, the electronic vana- 
dium spins slow down dramatically and thus induce an 
inhomogeneous static (on the NMR time scale) dipolar 
or transfered hyperfine field on neighboring ^^O sites and 
broaden the ^^O NMR line. Due to the fast decrease of 
signal intensity below Tn with decreasing temperature, 
detailed measurements of the temperature dependence of 
the line broadening in this temperature range were not 
pursued. 

The recovery of the ^^O longitudinal nuclear magne- 
tization AI (t) following the saturation pulses is a single 
exponential function at T > 100 K. This shows that all 
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Absorption spectrum of the O NMR signal for 
an " O-enriched powder sample of CaV2 04 at three differ- 
ent temperatures in an applied magnetic field H — 3.0 T. 
A strong inhomogeneous broadening is observed close to the 
magnetic transition temperature Tn = 78 K. The solid line 
at 77 K is a guide to the eye. 
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FIG. 5: Temperature T dependence of the ^'^O nuclear spin- 
lattice relaxation rate l/Tj [see Eq. ((T|] of ^^O-enriched pow- 
der CaV204 in applied magnetic fields H = 3.0 and 4.7 T. A 
strong peak is observed close to the antiferromagnetic transi- 
tion temperature Tn = 78 K, as labeled by the arrow. Inset: 
the stretching exponent P versus T. 



IV. 



NMR BELOW Tn IN SINGLE 
CRYSTALS OF CaV204 



the ^^O NMR transitions are saturated and thus there 
are no satellite peaks located away from the main peaks 
displayed in Fig. 31 consistent with a small quadrupole 
splitting as discussed above. Below 100 K, deviation from 
single exponential behavior is observed, indicating a dis- 
tribution of 1/Ti values. In order to extract a charac- 
teristic relaxation rate at all temperatures, we fitted the 
recovery curves to a stretched exponential function 



M{t) 
M{oo) 



= exp [-{t/T*f] 



(1) 



The relation of the parameters l/Tj* and f3 to the 1/Ti 
distribution underlying a stretched exponential function 
has been discussed in a recent paper ii2 Figure [5] shows 
the temperature dependence of l/T^ and P in H — 3.0 
and 4.7 T. At T > 100 K, the relaxation rate is almost 
temperature independent. Below 100 K, 1/Tj* exhibits a 
strong enhancement and reaches a peak at Tn = 78 K. 
Combining the above NMR results with magnetization 
studies of powder and single crystals,^ and with the mag- 
netic susceptibility data in Fig. [21 we identify the 1 /Tj* 
peak temperature as the temperature of an antiferromag- 
netic phase transition Tn. The peak in the nuclear spin- 
lattice relaxation rate at Tn results from an enhancement 
and slowing down of the electronic spin fluctuations at 
wave vectors close to the antiferromagnetic ordering vec- 
tor as the temperature approaches Tn from either side.^^ 



A. Spin Structure at 4.2 K 

In an external magnetic field H, the resonance fre- 
quency / of the ^^V nuclear spins is given by 



/ - |A(S> + /f|7v/27r, 



(2) 



where A is the hyperfine coupling constant between the 
nuclear spin and the vanadium electronic spins S, 7v is 
the gyromagnetic ratio of ^^V nuclear spins, and (S) de- 
notes the average electronic spin value in thermal equi- 
librium. In our experiments, the local field is much larger 
than the apphed field: H. Depending on 

whether the applied field component along (S) is paral- 
lel or antiparallel to (S), the resonance frequency shifts 
to higher or lower values, respectively. If CaV204 is a 
collinear antiferromagnet at low temperatures, where op- 
posite spin directions exist, upon application of an exter- 
nal field along the ordering axis, the zero-field spectrum 
should split into two peaks. These peaks should be sym- 
metrically displaced above and below the zero-field peak 
frequency. 

Figure [6] shows the field-swept spectra with the field 
at different angles in the a-b plane, measured with rf fre- 
quencies both higher and lower than the zero-field peak 
frequency /o = 236.7 MHz. All measurements in this 
subsection were performed on CaV204 crystal =f/=l. In 
contrast to a single peak above and a single peak below 
the zero-field peak frequency /o when an applied field is 
present, as expected for a collinear antiferromagnet, in- 
stead we see two peaks above /o and two peaks below 
/o in applied fields as shown by the vertical arrows in 
Fig. ini Whether each set of two peaks is resolved de- 
pends on the angle of the applied field in the a-b plane 
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FIG. 6: Field swept spectra with the applied magnetic field 
parallel to the a-b plane at rf frequencies of 231 MHz (left 
panels) and 243 MHz (right panels). The angles between the 
field and the average of the two projections of the two spin 
directions onto the a-b plane {S"^ in Fig. I10|l are labeled in 
each panel. The arrows indicate the positions of the peaks. 
S'm is approximately parallel to the crystallographic 6-axis. 
The spectra were measured at 4.2 K on crystal #1. 



as shown. We infer below that the spectra in Fig. [5] (and 
[8]) are consistent with a magnetic structure at 4.2 K that 
consists of two antiferromagnetic substructures, each of 
which is a coUinear antiferromagnetic arrangement where 
the angle between the ordering axes of the two substruc- 
tures is 19(1)°. The angle labeled in each panel of Fig. [6] 
is the angle between the applied field and the axis S'^, 
which is the average of the projections onto the a-b plane 
of the two spin ordering directions (see Fig. [TOl below). 

is approximately parallel to the 5-axis and is deter- 
mined by fitting the peak positions versus angle, as will 
be explained below. 

In our discussions of the ^^V NMR results, we assume 
that the applied magnetic field only shifts the NMR fre- 
quency without affecting the electronic spins. In fact, 
the ordered electronic moments can be tilted by the ap- 
plied field due to the presence of a torque. However, we 
can show that the tilting angle is indeed negligibly small. 
From the magnetization measurements,— at 4.2 K, the 
susceptibility x of single crystal CaV204 with applied 
field in the a, 6 or c directions is % ~ 0.003 cm^/mol, 



which corresponds to an induced moment of 0.005 /ie for 
each vanadium spin in a 2 T field. The tilting angle re- 
quired to produce such a moment is only 0.3°, assuming 
an ordered moment of 1.06 /ib-^ Both this angle and the 
induced moment are negligible to our studies. 

The spectra in Fig. [6] exhibit a two-peak structure 
when the field points away from the S'^ axis. When 
measured by field sweep at a fixed frequency /, each peak 
i {— 1,2) should appear at the field value 



peak i 



-/ocos ai 



9 • 2 

sm ai 



/(7v/2^), (3) 



where ai is the angle between the field and the respective 
electronic spin ordering direction of a magnetic substruc- 
ture noted above. The observed two-peak structures for 
/ > /o and / < /o strongly indicate the presence of 
these two different antiferromagnetic spin ordering direc- 
tions in the system. The ability to resolve the two peaks 
at the larger angles in Fig. [S] (and [5]) is related to the 
larger partial derivative of \dHpcaki/dai\ of Eq. ([3]) at 
the larger ai values for ai < Tr/2 rad. 

The two peaks in the same spectrum have different 
heights, as can be clearly seen in Fig.[S](b), (d), (g), and 
(i). The reason behind this difference is currently not un- 
derstood. It may be due to the change of the nuclear spin- 
spin relaxation times at different field values, since we fix 
the separation between the two rf pulses for echo gener- 
ation to be 18 lis. In fact, we observed a field-dependent 
oscillation of the spin echo intensity as a function of the 
separation between the two echo-generating pulses. Fig- 
ure [7] displays two nuclear spin-spin relaxation curves 
measured under identical conditions except for different 
external magnetic fields of 1.346 and 1.86 T, respectively. 
The oscillation pattern shows a clear field dependence. 
However, without knowing the detailed functional form 
of the decay curves, we cannot extrapolate the signal in- 
tensity back to zero pulse separation to correct for the 
nuclear spin-spin relaxation effect. 

We further note that the different heights of the two 
peaks cannot be attributed to the different percentage of 
spins in the two spin substructures. In such a scenario, 
the spins contributing to the left peaks in Figs, and 
(g) should contribute to the right peaks in Figs.IDJd) and 
(i) (see Fig. [TT]), and the left peaks should be higher in 
one orientation while lower in the other in Figs. [S]Jb) and 
(d), and in Figs.lHl^g) and (i), respectively. However, the 
spectra in Figs.IH^b) and (d), and in Figs.[6{g) and (i) are 
almost the same. The symmetry in the spectra with the 
field on opposite sides of S'^, such as in Figs, ^h) and 
(d) and in Figs. [Stg) and (i), indicates that the number 
of spins in the two substructures are the same. 

Since the above two-peak structure is observed with 
the field in the a-b plane, these measurements can only 
detect the difference of the spin projections of the two 
substructures onto the a-b plane. In order to determine 
whether or not the projections onto the b-c plane are also 
different, we also measured the spectra with the field in 
the b-c plane. Some representative spectra with the field 
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FIG. 7: Semilog plot of nuclear spin-spin relaxation curves 
at external magnetic fields of 1.346 and 1.86 T. The errors 
on the data points are smaller than the sizes of the symbols. 
The fields are parallel to the a-b plane and form an angle of 
31° from the S'm direction. Inset: the field swept spectrum 
measured at a fixed pulse separation of 18 /is in the same 
field direction. The vertical arrows show the field positions 
of 1.346 and 1.86 T. The measurements were performed at 
4.2 K on crystal #2. 



in different directions are displayed in Fig. [8l where the 
angles listed are described in the caption. As one can see, 
a two peak structure is still observed when the field is at 
a large angle from the S'^^ axis. However, the separations 
between the two peaks are smaller than in Fig.[6l indicat- 
ing a smaller angle between the two easy axis projections 
onto the b-c plane than onto the a-b plane. 

In order to study whether there exists canting and/or 
an imbalance in the number of spins in opposite direc- 
tions for each of the two ordered magnetic substructures, 
we compared the spectra with those measured with the 
field rotated by 180°. Figure [9] displays two spectra mea- 
sured at / = 222 MHz with the field parallel to the a-b 
plane and —31° and 149° away from the S'„^ direction, re- 
spectively. These two spectra are identical within exper- 
imental error, indicating the absence of spin canting and 
the same number of spins in opposite directions within 
each magnetic substructure. 

Thus we propose a model of the spin structure as shown 
in Fig. 1101 Various notations used in the model are ex- 
plained in the caption of Fig. 1101 There are equal num- 
bers of spins in the two antiferromagnetic substructures, 
each of which consists of coUinear antiparallel spins also 
with equal number. The plane defined by the two or- 
dering directions is parallel neither to the a-b nor the 
b-c plane. The average ordered moment direction Sm is 
approximately parallel to the 6-axis. This is consistent 
with single crystal anisotropic magnetization measure- 
ments versus temperature which showed that below Tn, 
the average easy axis of the magnetic structure is ap- 
proximately the 6-axisj^ Note that in our NMR study, 
we cannot determine the location in the lattice of the 
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FIG. 8: Field swept spectra with the applied magnetic field 
parallel to the h-c plane at rf frequencies of 231 MHz (left 
panels) and 243 MHz (right panels). The angles between the 
field and the average of the two projections of the two spin 
directions onto the h-c plane (S'^ in Fig. I10|l are labeled in 
each panel. The arrows indicate the positions of the peaks. 
The spectra were measured at 4.2 K on crystal #1. 
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FIG. 9: Comparison of two spectra measured at rf frequency 
/ = 222 MHz with and without the field direction reversed. 
The fields are parallel to the a-h plane and form angles of —31° 
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were measured at 4.2 K on crystal #1. 
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2 - 



FIG. 10: The proposed ordered spin structure in CaV204. 
There are two different antiferromagnetic ordering substruc- 
tures with equal numbers ol spins, each ol which has a 
collinear antilerromagnetic spin arrangement. AS, AS', and 
AS" are the angles between these two directions, and their 
projections on a-b and b-c planes, respectively. Sm , S'^ , and 
S!^ are the average of the two directions and their projections 
on a-b and b-c planes, respectively. 6'^ {O'h) and 6'^ {6'L) are 
the angles formed between a fixed arbitrary experimental ref- 
erence direction in the a-b {b-c) plane and the applied field H 
and S'jyi (S'/n), respectively. 



two different magnetic substructures. 

To extract the angle between the spin ordering direc- 
tions of the two substructures, we measured the depen- 
dence of the peak positions versus the field directions at 
rf frequencies of 231 and 243 MHz. The results are shown 
in Fig. [TT] For small angles between H and S'^ or S'''^, 
only one peak is observed. At larger angles, the positions 
of two peaks can be resolved. Two different symbols are 
used to represent the two different spin substructures, 
while for spectra with single peaks, a third symbol is 
used. Since the angle between the spin direction and the 
6-axis is much less than one radian, cosa^ in Eq. ([3]) for 
the field in the a-b plane can be approximated by 
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cosai,2 « cos(6»;, ± A6172 - O'h), 
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FIG. 11: Dependence of the peaks in the spectra at 4.2 K on 
the direction of the applied magnetic field, with the field in 
the b-c (top two panels) and a-b (bottom two panels) planes 
of crystal #1, where the rf frequencies are equal to 231 and 
243 MHz, respectively. For definitions of the angles 9'u, 9'^, 
O'/j, and O'm, see Fig. [TO] Circles and filled squares corre- 
spond to the two different spin ordering directions of the two 
magnetic substructures, respectively. The symbol * is used 
when the two peaks from the two spin directions overlap and 
only a single peak can be observed. The error in Hpcak is 
comparable to the size of the symbols unless shown explicitly. 
The solid and dotted lines represent the fits by the theoretical 
prediction in Eq. Q. 



where ai^2 are the angles between the field and the two 
spin directions Si and ^2, respectively. Similarly, for the 
field in the b-c plane, one has 

cos ai,2 ~ cos(C ± A6i'72 - 9'lj). (5) 

We fitted Eq. (O with cosai,2 given by Eqs. (g]) and ^ 
to the data in Fig. [TlJ The free parameters in the fit 
were /o, C C> : and AO". The best fit results are 
/o = 236.7(2) MHz, A0' = 18(1)°, and AO" = 6(1)°. 
The fits are shown in Fig. [TT] Since the angles between 
the spins and the 6-axis are much less than one radian. 



we have 

sin A9 w Vsin^ A6I' + sin^ AO", (6) 

from which one obtains the angle between the easy axes 
of the two magnetic substructures to be A9 = 19(1)°. 

In addition to the study of the angular dependence 
of the peak positions, we also measured their frequency 
dependences to further confirm the proposed spin struc- 
ture. Figures [T^ a) and (b) show the field swept spectra 
with the field H pointing along the S''^ direction, at rf 
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FIG. 12: (a) and (b): Field-swept ^^V NMR spectra at four 
different frequencies at 4.2 K. The frequencies are given under 
each spectrum in units of MIfz. The field is applied parallel 
to the S'm direction, (c): The frequency versus the peak field 
of the spectra. The solid lines are linear fits by Eq. (0 . The 
measurements were done on crystal #1 at 4.2 K. 



frequencies lower and higher than /o — 236.7 MHz, re- 
spectively. Note that when H \\ S'"j, the two magnetic 
substructures have the same peak positions of the spec- 
tra (see the zero-angle data in Fig. [TT|) . The peaks in 
Figs. [T^ a) and (b) both shift to higher fields when the 
frequency shifts further away from /o, respectively. The 
peak positions Hponk versus rf frequency are plotted in 
Fig. [iWcV The two sets of data points can be well fitted 
by the two linear equations 



/ = /o ± i/poak7/27r, 



(7) 



where /o is the peak frequency of the spectrum at zero 
applied field. A fit of Eq. ([7]) to the data gives '-f/2n — 
11.4(2) MHz/T, and fo = 236.7(1) MHz. Assuming 
Sm to be parallel to the 6-axis, the value of 7 should 
be 7/27r = (7v/27r) cos(A6l/2) = 11.07 MHz/T, where 
A9 = 19° is the above angle between the ordering direc- 
tions of the two magnetic substructures. This value of 
7/27r is very close to the above fitting result. 



B. Ordered Moment and Its Temperature 
Dependence 

In this subsection, we will study the temperature de- 
pendence of the vanadium ordered moment, which pro- 
vides evidence for an energy gap in the antiferromagnetic 
spin wave excitation spectrum, arising from anisotropy 
effects. Then we will discuss the value of the saturation 
vanadium spin moment at low temperatures. Measure- 
ments in this subsection were performed on crystal #2. 
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FIG. 13: Rf frequency / versus peak -ffpcak in field swept 
spectrum in crystal #2 at 4.2 K. The field is applied along 
the direction. The solid lines are fits with Eq. ([7]). 



The experiment was set up to allow field rotation in the 
a-b plane. By rotating the field in the a-b plane, the S'^^ 
direction (see Fig. [TU)) is identified as the direction along 
which the peak position in the spectrum is a minimum 
at a fixed rf frequency away from fo (see Fig. [TT|) . Af- 
ter identifying the S!^ direction, all subsequent measure- 
ments of the spectra were performed versus H at fixed rf 
frequencies with the field along the S!^ direction. 

With the field along S^, only a single peak is ob- 
served in the spectrum at each frequency (see Fig. I12p . 
In Fig. [131 we display the frequency dependence of the 
peak position at 4.2 K. Similar to the case of crys- 
tal #1 (Fig. [T^ . the data points can be well fitted by 
two straight lines. A fit of Eq. ([7]) to the data gives 
fo = 236.98(8) MHz and 7/27r = 11.3(1) MHz/T. This 
value of 7 is in agreement with the fitting value in crystal 
^1. However, the value of fo is slightly larger than in 
crystal #1. This slight difference may be due to sample- 
dependent differences. 

Figure [14] displays representative spectra measured at 
four different temperatures. For comparison between the 
different spectra, the a;-axis has been converted to the 
quantity / + Hj/2n, with 7/27r = 11.3 MHz/T. As 
the temperature increases, the signal intensity decreases 
rapidly and the spectra can only be measured below 45 K. 
In order to more accurately determine the peak position 
of the spectra, we fitted the original field swept spectra 
(with the X-axis being H) by a Gaussian function 

liH) =A + B CXp[-2(i/ - Hpeak)VM^'], (8) 

with A, B, iJpcak, and W as fitting parameters. The zero 
field peak frequency fo is then determined from 

/o = / ± i/pcak7/27r, 

where 7/27r = 11.3 MHz/T and the + and — signs cor- 
respond to the cases of / < /o and f > fo, respectively. 
In order to determine whether / < /o or / > /o, spec- 
tra were measured with at least two different frequencies 
at each temperature. With the correct choices of the -I- 
or — signs, the obtained fo values for different spectra 
as in Fig. [TJ] are the same within experimental error at 
each temperature. The final fo value is an average over 
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FIG. 14: Field swept spectra at different temperatures on 
crystal #2. The temperature and the rf frequency for each 
measurement are labeled in each panel. The solid lines are 
fits by Eq. ((Sjl to extract the peak positions. 



all calculated /o values for various spectra at the same 
given temperature. 

Figure [T2] shows the temperature dependence of /q. 
Since the temperature dependence of the hyperfine cou- 
pling constant can be ignoredjJ^ fo{T) is directly pro- 
portional to the ordered local moment. The ordered mo- 
ment is almost temperature independent at T < 15 K. 
In the spin wave theory of a three-dimensional antifer- 
romagnet without anisotropy, the initial decrease of the 
ordered moment with temperature T should follow a 
dependence. Fitting the data by a power law gives an 
exponent > 3.5, an unphysically large value (not shown). 
The temperature independence below T < 15 K thus in- 
dicates the presence of an anisotropy-induced energy gap 
for spin wave excitations^^ 

Before estimating the gap energy, we will first estimate 
the value of the saturation moment at T = K. The local 
moments reach their saturation value at T < 15 K. From 
the value of /o = 237 MHz at T < 15 K, one obtains a 
local field value of iJioc = /o/(7v/27r) = 21.2 T. In or- 
der to infer the value of local moment from the local field 
value, the hyperfine coupling constant A and the g-factor 
have to be determined. With known values of A and g, 
the ordered moment (/i^) is \{lJ-z)\ = ffMB-ffioc/^- The lo- 
cal field is dominated by the contact interaction through 
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FIG. 15: Temperature dependence of the NMR spec- 
tra peak position in zero applied field in CaV2 04 crystal 
#2. The dotted and solid curves are fits by Eqs. (|19p with 
one-dimensional spin wave dispersion and (|2U|I with three- 
dimensional dispersion, respectively. 



the polarized core electrons, which is approximately pro- 
portional to the number of unpaired electronic spins in 
the 3d orbitals with a proportionality constant of 12.5 T 
per unpaired electron to within 20%.^ However, beside 
the contact interaction, orbital effectsi^ and transfered 
hyperfine coupling with the neighboring ions^^ may 
also contribute significantly to the local field at the ^^V 
nuclear site. 

In the absence of a knowledge of the orbital effects 
and transfered interactions, we will estimate a possible 
range of the A/(g^B) value in CaV204 using the known 
values of A/{gfj.B) in other V'^+ compounds. The values 
of A/{gfiB) in LaVOg, YVO3, and V2O3 are 16.8, 20.8, 
and 15.5 T/fiB, respectively. For in AI2O3, EPR 
measurements gave A/{gfiB) — 13.3 T/hbi^ The range 
of A/{gfj,B) in the above four compounds is between 13.3 
and 20.8 T/^b- Assuming A/{gfiB) in CaV204 hes in 
the same range, the low temperature ordered vanadium 
moment in CaV204 is thus in the range of 1.02-1.59 fiB- 
This \{fiz)\ range is consistent with the value 1.06(6) /^b 
extracted from the previous neutron powder diffraction 
study ii 

Next we estimate the energy gap for the spin wave ex- 
citations. From the neutron diffraction studies,^ we know 
that the spins reverse their ordering direction upon mov- 
ing along the c-axis. Therefore the spin structure in the 
ordered state within a zig-zag chain should look as in 
Fig. 1161 where we ignore the possible small misalignment 
of 19° between spins in the two magnetic substructures 
discussed above. Because each spin in one leg of the 
chain couples by the same exchange constant Ji to two 
spins in the other leg that are ordered in opposite di- 
rections, we expect the effective coupling between the 
two legs within a zig-zag chain is greatly reduced in the 
ordered state. As a result, we will consider the nearest- 
neighbor interactions within the leg J2 as the only domi- 
nant magnetic interaction and treat the effect of interleg 
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FIG. 16: Zig-zag spin structure in CaV204. Due to the alter- 
nation of the spin directions along the c-aixs, the interaction 
between the spins in the two legs of the zig-zag chain are 
essentially decoupled. The possible misalignment of 19° be- 
tween spins in the two legs of the zig-zag chain is ignored. Ji 
and J2 are the nearest-neighbor inter and intra-leg exchange 
interactions, respectively. 



interaction within a zig-zag chain as a weak interchain 
interaction. For simpUcity, we will use a single exchange 
constant J' to characterize the effect of the interchain in- 
teractions. To include the effect of anisotropy, we assume 
a single ion anisotropy in the system with a single direc- 
tion of easy axis. Then, the Hamiltonian can be written 
as 



2 ^ J' Si ■ Sj+i, 



(9) 

where K is the anisotropy constant, the index i runs 
through the spins in one leg of the chain, and the summa- 
tion (i, j) runs through all interleg and interchain nearest- 
neighbor pairs. 

The Fourier transform of the exchange interactions is 

J{q) = J{rij)cxp{-iq ■ r,j), 



where ry connects two spins in opposite sublattices, 
J{rij) > is the nearest-neighbor exchange constant, 
and the index j runs through all the nearest-neighbor 
spins of spin i in the opposite sublattice (each sublattice 
consists of spins in the same direction). The spin wave 
dispersion relation is given by^S 



Eg ^ {[2SJ{0) + KS]^ - [2SJ{q)]^} 



1/2 



where we ignored interactions between spins in the same 
sublattice. The spin wave gap value is given by the value 
of Eg at q = 0. In the limit of small anisotropy K <^ 
2^7(0) w 4J2, the gap energy is given by 



Eg - 2S[JiO)K]^^\ 



(11 



In the spin wave theory, the decrease of sublattice 
magnetization is due to the thermal activation of spin 
wave excitations. In the above bipartite antiferromag- 
netic system 



20 



V 



2J{Q)S + KS 



(12) 



where the integral is limited to the first Brillouin zone 
of one sublattice, V is the sample volume per sublattice 
site, and 

, > 1 



is the number of thermally excited antiferromagnetic 
magnons. 

The dispersion relation in Eq. pUj) depends on the spin 
lattice structure and the exchange interactions J{r). For 
a quasi one-dimensional chain with interchain coupling 
J', at temperatures T ^ J' /ks, one can ignore the dis- 
persion perpendicular to the chain direction. Then for 
small values of qc, which is the q vector component along 
the chain, one can perform a Taylor series expansion of 
1^(9)1' as 



(13) 



where I is the nearest-neighbor distance within the leg. 

At T <C 0/2/^8, only spin waves at small qc values have 
significant contributions to the integral in Eq. (|12p , so one 
can change the limits of integral for qc in Eq. (|12p to ±00. 
The small q approximation is valid only at temperatures 
where 1 - |(S'^(r))|/|(S'^(0))| < 0.1^ which is satisfied 
within our experimental temperature range. Substituting 
Eq. (|13p into Eq. (llOp . and changing the limits of integral 
for qc in Eq. (fT^ to ±00, one obtains in the limit of small 
anisotropy K < 2J{0) and T <^ Eq 



1 - 



/2 



where 



B : 



2tt 



0.80. 



(14) 



(15) 



Equation is valid at temperatures J'/Zcb <C T <C 
J2/fcB. In CaV204, J' /kB might fall within the experi- 
mental temperature range in Fig. [15] (1.5 < T < 45 K). 
So it is useful to consider the other limit of T ^ J' /k-Q <C 
J2/fcB, where a three dimensional dispersion is more ap- 
propriate. Applying a small q approximation, one obtains 

(10) \j{q)\' « 1^(0)1^ {1 - ^\v/2f'' + g?)] } , 

(16) 

where for simplicity, we assumed an isotropic dispersion 
in the a-b plane, i] is a geometrical factor of order one 
which depends on the spin structure^-^ and j' is of the 
order of J'/J2- By combining Eqs. and (fTC)) . 

and changing the three limits of integrations to ±00 in 
Eq. (fT^ . then instead of Eq. HH), we have in the limit 
of small anisotropy K < 2 J^(0) and T < Eq^ 



(SAT)) 



Be-^^l^^'^{EGlkBT)- 



■3/2 



where 



B 



V2a 



(17) 



(18) 
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witha = K/[2J{Q)]. 

Using Eqs. ^1]) and pT|) and the relation 

/o(r)//o(o) = |(5,(r))|/|(5,(o))|, 

one obtains the variation of the zero-field '''^V NMR res- 
onance frequency /q as (for ID) 



k{T) - /o(0) \l - Be-^<^"'-^{EG/k^T)-^ 



/2 



(19) 



or (for 3D) 

k{T) - /o(0) \l - Be-^«/'=-^(£;G/fcBr)-3/2l , (20) 



depending on whether a one-dimensional (ID) or three- 
dimensional (3D) dispersion is used for Eq. We fitted 
Eqs. in]) and (gni) to the /o(r) versus T data in Fig.IH] 
at T < 45 K with /o(0), B, and Eq as free parame- 
ters. The best fit results are /o(0) = 237.04(5) MHz, 
Eg = 98(5) K, B = 0.51(7) for the ID dispersion with 
Eq. dm), and /o(0) = 237.08(6) MHz, Eg = 64(5) K, and 
B = 0.27(6) for the 3D dispersion with Eq. The 
best fit curves are shown as the dotted and solid curves 
in Fig. [iSl respectively. Since Eqs. (fT9|) and ([20|) are 
derived under the two limiting conditions of T 3> J' /fee 
and T <C J'/fce, respectively, one can expect that the ac- 
tual Eg value might lie somewhere between 64 and 98 K. 
Given a value of Eg, we can make a rough estimate of the 
anisotropy constant K. From the magnetization study, 
one estimates that the intrachain nearest-neighbor ex- 
change constant to be of the order of J2/A;b ~ 200 Ki^ 
Taking S ^1,Eg = ?>1 K, and J{Q) w 2J2, we thus have 
from Eq. ^ that K/k^ - 4 K. 

The above fitting value of i? = 0.51(7) from ID dis- 
persion is similar to the calculated value of 0.80. The 
fitting value of i? = 0.27(6) from the 3D dispersion 
constrains the value of J' in the 3D model. Taking 
a « K/AJ2 ~ 0.01, 7? - 1, and B = 0.27, then from 
Eq. (dll) one has / - 0.01 and J'/fce j'J2/kB ^ 2 K. 
However, this value of J' seems inconsistent with the 



initial assumption of T <C J'/^b required for the three- 
dimensional model to be valid. Therefore, Eq. of the 
ID model might provide a better approximation to the 
/o(r) data. 



V. SUMMARY AND CONCLUSIONS 

We have presented ^^O and ^^V NMR results on the 
zig-zag spin chain compound CaV2 04. The strong inho- 
mogeneous broadening and a peak in the nuclear spin- 
lattice relaxation rate versus temperature of ^"^0 NMR 
confirm the presence of an antiferromagnetic phase tran- 
sition at 78 K in a powder sample. The crystals we stud- 
ied have Tn = 69 K. ^ly NMR in the ordered state 
of crystals reveals the presence of two antiferromagnetic 
substructures at 4.2 K, each of which is coUinear and 
which form an angle of 19(1)° between them with the av- 
erage direction approximately parallel to the &-axis. The 
origin and location in the lattice of the different spin sub- 
structures remain unknown. However, we speculate that 
the two magnetic substructures are associated with the 
two inequivalent V'^'*' S = 1 zig-zag spin chains in the 
orthorhombic crystal structure, respectively. Magnetic 
neutron diffraction studies are required to further char- 
acterize the magnetic structure below Tn. The temper- 
ature dependence of the zero-field resonance frequency 
at low temperatures suggests the presence of magnetic 
anisotropy and an energy gap in the spin wave excitation 
spectrum. The energy gap is estimated from spin wave 
theory to be between 64 and 98 K. 
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